Halo Star Lithium Depletion 



M.H. Pinsonneault 1 , T.P. Walker 1)2 , G. Steigman 1>2 , and Vijay K. Narayanan 1 



Received ; accepted 



00 



> 

m 
o 

m 

o 

oo 

On 

S3 

9 L,: 
6 

CO 

C3 



X 



1 Department of Astronomy, The Ohio State Univerity 
2 Department of Physics, The Ohio State University 



- 2 - 



ABSTRACT 

The depletion of lithium during the pre-main sequence and main sequence phases of stellar 
evolution plays a crucial role in the comparison of the predictions of big bang nucleosynthesis 
with the abundances observed in halo stars. Previous work has indicated a wide range of possible 
depletion factors, ranging from minimal in standard (non-rotating) stellar models to as much as 
an order of magnitude in models which include rotational mixing. Recent progress in the study 
of the angular momentum evolution of low mass stars (Krishnamurthi et al. 1997a) permits the 
construction of theoretical models which reproduce the angular momentum evolution of low 
mass open cluster stars. The distribution of initial angular momenta can be inferred from stellar 
rotation data in young open clusters. In this paper we report on the application of these models 
to the study of lithium depletion in main sequence halo stars. 

A range of initial angular momenta produces a range of lithium depletion factors on the 
main sequence. Using the distribution of initial conditions inferred from young open clusters 
leads to a well-defined halo lithium plateau with modest scatter and a small population of 
outliers. The mass dependent angular momentum loss law inferred from open cluster studies 
produces a nearly flat plateau, unlike previous models which exhibited a downwards curvature 
for hotter temperatures in the 7 Li - T off plane. The overall depletion factor for the plateau 
stars is sensitive primarily to the solar initial angular momentum used in the calibration for the 
mixing diffusion coefficients. 

The 6 Li/ 7 Li depletion ratio is also examined. We find that the dispersion in the plateau 
and the 6 Li/ 7 Li depletion ratio scale with the absolute 7 Li depiction in the plateau and we use 
observational data to set bounds on the 7 Li depletion in main sequence halo stars. A maximum 
of 0.4 dex depletion is set by the observed dispersion and 6 Li/ 7 Li depletion ratio and a minimum 
of 0.2 dex depiction is required by both the presence of highly overdepleted halo stars and 
consistency with the solar and open cluster 7 Li data. The cosmological implications of these 
bounds on the primordial abundance of 7 Li are discussed. 
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Introduction 



The status of big bang nucleosynthesis (BBN) as a cornerstone of the hot big bang cosmology rests 
on the agreement between the theoretical predictions and the primordial abundances of the light elements 



deuterium (D), helium-3 ( 3 He), helium-4 ( 4 He), and lithium-7 ( 7 Li) inferred from observational data (Yang 
et al. 1984| ; |Walker et al. 1991 ; |Copi et al. 1995a ). Comparisons of this type often rely on models of galactic 



chemical and/or stellar evolution in order to associate the observed abundances of these light elements at 
the present epoch with their predicted primordial values. Recently the confrontation between prediction 
and observation has come under stricter scrutiny as it appeared that the primordial abundance of deuterium 
inferred from ISM/Solar data was smaller than its predicted abundance - the value of which follows from 
requiring that the standard big bang model predictions for 4 He are in good agreement with the abundance 
inferred from observations of metal-poor extragalactic H n regions ( [Hata et al. 1995| ; |Copi et al. 1995b| ; Hatd 



et al. 1997). Very recent measurements of the deuterium abundance along the lines-of-sight to high-redshift 
QSOs( Buries fc Tytler 1997 ), along with a reanalysis of the 4 He data in light of new observations (Olive 



Skillman, & Steigman 1997), increase the tension between prediction and observation. 



The role of Li in these comparisons has been minor due to the large uncertainty in the estimate of 
the amount of lithium destruction during the lifetimes of Popn (i.e., metal poor) halo stars which could 
accomodate primordial lithium abundances ranging from the observed plateau value 3 (no depletion) up to 
a factor of ten larger. This range in possible depletion factors is equally compatible with primordial lithium 
abundances which correspond to either "low deuterium" (which favors a primordial lithium abundance a 
factor of 3 higher than the plateau value), or the observed 4 He (which favors the plateau value). 

The arguments in favor of minimal lithium depletion are the flatness of the Popn lithium abundance 
plateau (at low metallicity and high temperature) with respect to both metallicity and temperature and 
the low dispersion in the lithium abundance at fixed T c ff. This was generally consistent with "standard" 
stellar models (i.e., models without rotation which burn lithium via convection during pre-main sequence 



evolution) which predicted little depletion (Dcliyannis et al. 1989). There are however some specific areas 
of disagreement in the comparison of the halo star data with standard models. The observed dispersion 



may be greater than that predicted (Deliyannis, Pinsonneault, & Duncan 1993; Thorburn 1994; Ryan ei 



3 [ Thc lithium abundances of Popn stars with T e g > 5800K and [Fe/H] < —1.3 is nearly independent of 
metallicity ([Li] = 2.25 ± 0.1, where [X] = 12 + \ogyx with yx the number ratio of X to hydrogen) and, 
hence, is referred to as a "plateau" . 
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al. 1996) or not (Molaro, Primas, & Bonifacio 1995; Spite et al. 1996). There is a small population of 
highly overdepleted stars which appear normal except for 7 Li (Thorburn 1994; Norris et al. 1997j ), and 
the trends with metal abundance appear to conflict with expectations from the models (Thorburn 1994, 



Chaboyer fc Demarque 1994 ). However, the overall agreement between standard stellar evolution theory 
and observations in halo stars is good enough that nonstandard models would probably not be invoked 
to explain this data set. Therefore many investigators have adopted the reasonable assumption that the 
observed Popn lithium abundances are close to the primordial value. 

The overall properties of 7 Li depletion in standard models have been extensively studied (for a review 



see Pinsonncault 1997 ) and there are some model independent predictions: There is partial 7 Li depletion 
in the pre-main sequence (pre-MS), little or no main sequence (MS) depletion for stars hotter than about 
5500 K, and there should be little or no dispersion in the 7 Li abundance at a given T e ff in clusters. 
Unfortunately, the observational data obtained from halo stars does not serve as a good diagnostic for 
these models since the initial abundance is not known and we have no information on the history of the 
7 Li abundance. Instead, one can look to the Sun and open clusters, systems with a nearly uniform initial 
lithium abundance ([Li] = 3.2 to 3.4) where detailed abundance data is available as a function of mass, age, 
and composition. These Popi data provide stringent tests of theoretical models and it has been known for 



quite some time (e.g., see Weyman fc Sears 1965 and Zappala 1972 ) that standard models fail these tests. 
The disagreement between the data and standard models has increased as both the observational data and 
the standard theoretical models have improved. The open cluster data provide clear evidence that the 7 Li 
abundance decreases with increasing age on the MS, contrary to the standard model prediction that the 
convective zones of these hot (> 5500K) stars are not deep enough to destroy 7 Li on the MS. In addition, 
there is strong evidence for a dispersion in abundance at fixed T e g and unexpected mass-dependent effects, 
such as the strong depletion seen in mid-F stars relative to both hotter and cooler stars (the "Li dip" of 



Bocsgaard fc Trippico 1986 - for extensive reviews of this issue, see Pilachowski fc Hobbs 1988 , Michaud & 



Charbonneau 199l| , poderblom 1995[ p3alachandran 1995[ and [Pinsonncault 1997| ). For standard models, 
these features cannot be explained by variations in the input physics (e.g., opacities); rather, they indicate 
the operation of physical processes not usually included in standard models. By extension, these processes 
could also be operating in the Popn stars. 

Another potentially useful diagnostic of lithium depletion can be found in globular clusters. Clusters 
provide samples which are homogeneous in age and composition, so one would expect a smaller dispersion 
in a cluster sample than in a field star sample. This is certainly the case for Popi stars: Popi field stars do 



- 5 - 



show a larger dispersion than open cluster stars (Lambert, Heath, & Edvardsson 1991; Favata, Micela, & 



Sciortino 1996). However, recent observations of Li in globular cluster subgiant and turnoff stars show the 
opposite trend and create some serious difficulties for the minimal 7 Li depletion scenario for metal poor 



stars ( Deliyannis, Bocsgaard, fc King 1995 , Boesgaard et al. 1997 ; Thorburn et al. 1997 ). There are clear 
star-to-star differences in M92; out of 3 stars observed with very similar color, 2 have lithium abundances 
well below the plateau value. In NGC 6397, 20 stars near the turnoff were observed. For 7 with identical 
B-V color, there is a scatter of a factor of 2-3 in lithium abundance. Therefore a standard model treatment 
of the halo field stars must explain why the lithium abundances in globular clusters, but not in field stars, 
are anomalous. 

The disagreement between the Pop i lithium abundances and the predictions of standard stellar models 
has stimulated investigation of nonstandard stellar models. The most prominent explanations are mixing 



(either from rotation or waves), microscopic diffusion, and mass loss (Pinsonneault 1997; Michaud & 



Charbonneau 1991). For Popi G and K stars, microscopic diffusion and mass loss are not likely explanations 



(Michaud & Charbonneau 1991, Swenson & Faulkner 1992), which leaves mixing as a logical candidate. 
Mixing induced by rotation can explain many of the overall properties of the Pop i data. Rotational mixing 
can, on relatively long timescales, mix material from the base of the convective zone to interior regions 
of the star where lithium can be burned; both the rotation rate and the 7 Li depletion rate decrease with 
increased age; and a distribution of initial rotation rates can produce a distribution of 7 Li depletions and 
thus a dispersion in abundance at fixed T e g. 

Although models with rotational mixing have the qualitative properties needed to explain the Pop 
I and Pop II data, two classes of objections have been raised: 1) the uniqueness of the solutions and 
adequacy of the physical model; and 2) discrepancies in the quantitative comparison of observation and 
theory. Rotational models require an understanding of the angular momentum as a function of radius 
and as a function of time. There have been persistent difficulties in reproducing the surface rotation 



rates as a function of mass and time in open clusters (Chaboyer et al. 1995a, b; Keppens, MacGregor, & 



Charbonneau 1995; Bouvier 1995). In addition, models with internal angular momentum transport from 
hydrodynamic mechanisms, such as the ones used in this paper, predict more rapid core rotation in the Sun 



than is compatible with helioscismology data ( Chaboyer et al. 1995a , Krishnamurthi et al. 1997a (KPBS), 
Tomczyk, Schou, fc Thompson 1995| ). Either of these difficulties could potentially affect the degree of 
mixing in the models. 
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Models with rotational mixing predict that a range of initial angular momenta will generate a range of 
lithium depletion rates and therefore a dispersion in abundance among stars of the same mass, composition, 
and age. However, the difficulty in reproducing the angular momentum evolution of low mass stars has made 
quantifying the predicted dispersion difficult, and therefore only qualitative estimates have been made (e.g., 
see Pinsonneault, Kawaler, fc Demarque 1990| (PKD), |Chaboyer fe Dcmarquc 1994 , |Chaboyer, Demarque, fc 
Pinsonneault 1995b, and pharbonnel, Vauclair, fc Zahn (1992) ). There have also been mismatches between 



the mean trend inferred from the data and from all classes of theoretical models. Chaboyer & Demarque 
(1994), for example, concluded that the observed mean trend of lithium abundance with T e g in metal-poor 
halo stars was not reproduced in any class of models, be they standard, rotational, or with diffusion. There 
are also cases, such as the cool stars in young open clusters, where the observed dispersion is large but the 
rotational models predict little, if any, dispersion. 

The case for or against standard model lithium depletion is not as clear cut when we consider halo stars. 
It is argued that models with significant rotationally induced depletion could not produce a flat plateau 



with limited dispersion about the mean plateau abundance (e.g., Bonifacio & Molaro 1997). Furthermore, 
it was expected that such models would destroy far too much 6 Li ( [gtcigman et al. 1993 ) in conflict with the 
observation of 6 Li in HD 84937 (see section 4.2). On the other hand, standard (convective burning) models 
show lithium depletion trends contrary to the Popi data. At minimum the mechanisms responsible for the 
Pop i 7 Li pattern need to be identified and shown to not affect Popn stars. It is difficult to construct a 
model where the nonstandard effects completely cancel for Pop n stars while still retaining consistency with 
the Pop i data. Models with microscopic diffusion predict modest 7 Li depletion in plateau halo stars. In 
general these same models predict that the timescale for changes in the surface Li abundance decreases 
with increased mass; this produces a downward curvature in the 7 Li-T e gf relationship at high temperatures 



which is not observed. Mass loss can counteract this trend ( Swenson 1995 ; Vauclair fc Charbonncl 1995 ). 
In either case, the net effect is modest depletion at the 0.2 dex level. Rotational mixing can also suppress 



diffusion (Chaboyer & Demarque 1994, |Vauclair fc Charbonncl 1995[) . The cancellation of different effects 
still requires at least some mean 7 Li depletion in the halo stars. In a recent preprint Vauclair & Charbonncl 



(1998) note that although the surface 7 Li abundance varies strongly with mass in diffusive models, the 



peak subsurface abundance does not (note that the height of the peak is not preserved in models which 
include mixing). They then use the height of that peak to constrain the absolute 7 Li abundance, arguing 
that the appropriate mass loss strips each star down to the region that contains this peak 7 Li abundance. 
The survival of 6 Li is also problematic in unmixed models with sufficient mass loss to expose the peak 
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abundancc; in a model of ours which can be compared to Figure 1 in Vauclair & Charbonncl (1998) 



'Li 



drops to half its surface value in the outer 0.01 solar masses, well below the comparable mass content in the 



7 Li preserving region of 0.02 solar masses and also the peak in the Vauclair fc Charbonncl (1998) model. 
Further, it is disturbing to use a class of models which do not include a fundamental characteristic of stars, 
namely, rotation, when that property has been shown to be capable of affecting the issue being studied. 

In our view the most serious source of uncertainty in models with rotational mixing has been the 
understanding of the angular momentum evolution of low mass stars. We will show that many of the 
difficulties in reconciling observation and theory are resolved when models that are consistent with the 
rotation data in open clusters are used. In order to make definitive predictions of the amount of rotational 
mixing, it is necessary to follow the stellar angular momentum histories which requires knowledge of 
the initial distribution of rotation velocities along with an angular momentum loss law. There has been 
significant recent progress in constructing theoretical models consistent with the rotation data in young 



open clusters (KPBS; Collier-Cameron & Li 1994; Keppens, MacGrcgor, & Charbonneau 1995; Bouvier, 



Forestini, & Allain 1997; Allain 1997). With the latest generation of models, we can both infer the 



distribution of initial conditions and place strong constraints on the surface rotation as a function of time. 
This enables us to quantify the expected dispersion in lithium abundance and greatly reduces the sensitivity 
of the model predictions to uncertainties in the input physics. Those models which accurately reproduce 
the lithium abundances observed in open clusters can then be used to predict lithium depletion in halo 
stars. The general trend is that these "open cluster normalized" rotation models for halo stars predict more 
lithium depletion than do the standard models but less depletion than that predicted by earlier studies of 
rotational mixing with a less sophisticated treatment of angular momentum evolution. 

Our goal in this paper is to constrain the amount of lithium depletion in Pop n stars using several 
observables: (1) an estimate of the halo star initial rotation rate distribution as derived from open clusters, 
(2) the absence of large dispersion in the observed lithium abundances of the Popn "plateau" stars, (3) 
the 6 Li abundance and/or the 6 Li/ 7 Li abundance ratio in HD 84937. We argue that, individually and in 
combination, these observables point to Popn halo star lithium depletion of at least 0.2 dex (following 
from the open cluster data) but no more than 0.4 dex (a consequence of the narrowness of the plateau 
and of the 6 Li considerations). Based on the observational data we can then infer the primordial lithium 
abundance and compare and contrast it with that predicted by standard BBN for consistency with the 
inferred primordial abundances of D and/or 4 He. 
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2. Method and Comparison with Previous Models 
2.1. Angular Momentum Evolution and Rotational Mixing 

In standard stellar models, lithium depletion is a strong function of mass and composition; it also 
depends on the input physics, particularly the opacities and model atmospheres used to specify the 
surface boundary condition. In our models the standard model physics is the same as described in KPBS 



( Krishnamurthi et al. 1997a ) : namely, interior opacities from OPAL ( Roger fc Iglesias 1992 ) , low-T opacities 
and model atmospheres at solar abundance from Kurucz (1991a ^); nuclear reaction rates from Bahcall, 



Pinsonneault fc Wasserberg (1995)] , including the Caughlin & Fowler (1988) 7 Li(p, a) cross-section; Yale 
EOS; a solar calibrated mixing length, and Y ~ 0.235. Chaboyer & Demarque (1994) noted that unusual 
Li-T e ff trends for cool metal-poor stars occurred in models using the Kurucz atmospheres and could be 
lessened by using a grey atmosphere. We therefore ran halo star models with grey atmospheres and a 
mixing length of 1.25 as obtained for solar models constructed under the same assumptions. 

Stellar models which include rotational mixing require additional input physics beyond standard stellar 
models. The important new ingredients include: 



1. a distribution of initial angular momenta 

2. a prescription for angular momentum loss 

3. a prescription for the internal transport of angular momentum and the associated mixing in radiative 
regions 

4. the impact of rotation on the structure of the model. 



There have been important changes in the treatment of the first three ingredients since the study of 



rotational mixing in halo dwarfs by [Pinsonneault, Deliyannis, fc Demarque (1992; PDD)| ; in this section we 
discuss the ingredients of the current set of models and compare them with previous work. The treatment 
of angular momentum evolution in this set of models is the same as KPBS. The structural effects of rotation 



have been computed using the method of Kippenhahn & Thomas (1970), and are small for low mass stars 
which experience angular momentum loss. 
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2.1.1. Initial Angular Momentum 



The studies of lithium depletion in the presence of rotationally induced mixing by Pinsonneault ei 



al. (1989; PKSD) (the Sun), PKD (open cluster stars), and PDD (halo stars) all used a range of initial 
rotation rates in the pre-MS from the measured rotation velocities of young T Tauri stars (10-60 km/s) at 
a typical reference age of 1 Myr. No interaction between accretion disks and the protostar was included; 
the entire range of initial conditions was generated from the range in the initial rotation velocity. This 
caused some difficulty in reproducing the observed range in rotation (a factor of 20 in young open clusters), 
especially since the higher angular momentum loss rate in rapid rotators acted to reduce the range in initial 
rotation rates. This approach was also used by Chaboyer et al. (1995a,b) for the Sun and open cluster stars 
respectively and by Chaboyer & Demarque (1994) for halo stars, although these papers used a different 
angular momentum loss law than the earlier work. 

As protostars evolve their moment of inertia decreases; this would produce an increase in rotation rate 
with decreased luminosity. In contrast, the rotation periods of classical T Tauri stars appear to be nearly 
uniform - they do not scale with luminosity as they would if the stars were conserving angular momentum as 



they contracted towards the main sequence (Bouvier et al. 1993, Edwards et al. 1993). Weak- lined T Tauri 



stars, which are thought to be protostars without significant accretion disks, have much shorter rotation 
periods. This behavior would be expected if there were an accretion disk regulating the rotation of the 
central object and spinup only occurred when the disk is no longer magnetically coupled to the protostar 



( Konigl 1991 ; Cameron, Campbell fc Quaintrcle 1995 ; Kcppens, MacGregor, fc Charbonneau 1995 ). In this 



revised picture there is a constant surface rotation rate while there is a sufficiently massive disk around the 
central object; stars which detach from their disks earlier experience a larger change in angular velocity as 
they contract to the main sequence than those which detach from their disk later (and are only free to begin 
to spin up when their moment of inertia is smaller) . This disk- locking hypothesis leads to a larger range in 
rotation rates for models of young open cluster stars, in better agreement with the data, and explains the 



very slow rotation of the majority of young stars (Bouvier, Forestini, & Allain 1997; Cameron, Campbell 



fc Quaintrcle 1995| ; Keppens, MacGregor, fc Charbonneau 1995 ; KPBS). These differences in the initial 



conditions lead to interesting consequences for the lithium depletion pattern. Although the initial rotation 
periods are comparable to those used in the earlier studies, the inclusion of star-disk coupling implies that 
the MS angular momentum is much lower for models with long disk lifetimes than it would be if the disks 
were absent. Disk lifetimes of 3 Myr are needed to match the rotation rates of the slow rotators in young 
clusters. The mean depletion for the majority of stars is therefore significantly lower than in previous 
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studics because most young stars are slow rotators and the degree of mixing decreases for models with lower 
MS angular momentum. 



2.1.2. Angular momentum loss. 

PKSD, PKD, and PDD all used a loss law of the form dJ/dt oc oj 3 , where oj is the angular velocity 
( Kawaler 198§| ). This implies greater angular momentum loss for rapid rotators than for slow rotators. 
However, models of this type do not reproduce the rapid rotators seen in young open clusters; the spin down 
of high angular momentum objects is too severe (PKD). On both theoretical and observational grounds, a 
more realistic functional form for the loss law is 



and 



dj n . 

— OCCJ (W < Wcrit), 

dJ 9 



dt 



OC 0JUJ r 



(1) 



(2) 



( [MacGregor fc Brenner 199lfc |Cameron, Campbell fc Quaintrelc 1995| ; |Mestel 1984] |Barnes fe Sofia 1996| ). 
Here w C rit is the angular velocity at which the angular momentum loss rate saturates; w cr it can be estimated 
from several observable quantities, thought to be correlated with surface magnetic field strength, which 



saturate at 5-20 times the solar rotation period (see Patton fc Simon 1996, Krishnamurthi et al. 1997b 



for reviews and discussions of recent work). A loss law which saturates at high rotation rates allows rapid 
rotation to survive into the early main sequence phase; a loss law of this form was adopted by Chaboyer et 
al. (1995a, b) and Chaboyer & Demarque (1994). 

The observational data indicate that the duration of the rapid rotator phase is a function of mass 
in the sense that rapid rotation survives longer in lower mass stars. There is also indirect observational 



evidence for a mass-dependent saturation threshold (Patton & Simon 1996; Krishnamurthi et al. 1997b). 
Models where the saturation threshold scales with the convective overturn time scale can reproduce the 
observed mass dependent spin down pattern. We adopt the mass-dependent saturation threshold of KPBS, 



T conv(©) /„x 

Wcrit = Wcrit(O) 7T, (3) 

7"conv(*J 

where the convective overturn time scale r conv as a function of ZAMS T c g was taken from the 200 Myr 



isochrones of Kim & Demarque (1996; KD). KD only considered models at solar abundance, and we would 
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need to recalibrate the KPBS models if we used a time and abundance dependent overturn timescale for the 
models. A comparison of the convection zone depth (as a function of T e fi, Z, and age) in the tables of PKD 
and PDD indicates that at young ages the convection zone depth at a given T e ff , and therefore the overturn 
time scale, depends only weakly on metal abundance. The value of w cr j t is important primarily in the early 
main sequence. We therefore evolved standard halo star models to an age of 200 Myr and used the same 
Tconv as a function of T e ff for the metal-poor models as reported by KD for the solar abundance models. 
We will consider models with a mass, composition, and time dependent overturn time scale in a future 
paper in preparation ( Narayanan et al. 1998| ). A loss law with a mass-dependent saturation threshold acts 
preferentially to suppress rapid rotation in the more massive stars, unlike the constant saturation threshold 
adopted by Chaboyer & Demarque (1994). This reduces both the absolute depletion and the dispersion of 
lithium abundances in hotter halo plateau stars relative to cooler plateau stars. 



2.1.3. Angular momentum transport. 



The treatment of angular momentum transport is the same as in KPBS; a detailed review of the time 
scale estimates can be found in Pinsonneault (1997). We consider internal angular momentum transport 
by hydrodynamic mechanisms alone, and do not include potentially important mechanisms such as gravity 
waves (Kumar & Quartaret 1997; Zahn, Talon, & Matias 1997) and magnetic fields QCharbonneau & 



MacGregor 1992, 1993; Barnes, Charbonneau, & MacGregor (1998)). In previous work the degree of 
mixing was found to be insensitive to the assumptions about internal angular momentum transport (PKD; 



PDD; Chaboyer et al. 1995a , b). Charbonnel, Vauclair, fc Zahn (1992) and Richard et al. (1996) obtained 



similar depletion patterns in models with solid body rotation, and Barnes, Charbonneau, fc MacGregor 



(1998) reported similar results in models with magnetic fields and hydrodynamic mechanisms. This can 



be understood because the diffusion coefficients for angular momentum transport are actually largely 
determined by the balance between the surface boundary condition (the applied torque) and the flux 
of angular momentum from below a given shell. Theory gives an estimate of the ratio of the diffusion 
coefficients for mixing to those for angular momentum, which is related to the existence of anisotropic 



turbulence in stellar interiors ( Chaboyer and Zahn 1992 ). The diffusion coefficients for mixing are calibrated 
on the Sun, but our lack of information on the rotational history of the Sun requires an exploration of 
models with varying solar initial conditions. 
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3. Results 

3.1. Initial Conditions and Angular Momentum Evolution 

We have taken the Pleiades vsmi data used in KPBS and applied a statistical correction of 4/ir to 
all of the data points to correct for inclination angle effects. A grid of solar composition models with a 
range of accretion disk lifetimes was constructed and evolved to the age of the Pleiades. We interpolated 
in the grid to map the distribution of measured rotation velocities onto a distribution of accretion disk 
lifetimes. A histogram of the number of stars as a function of disk lifetimes is presented in Figure la. This 
distribution of initial conditions was used to infer the distribution of lithium depletion factors for different 
solar disk lifetimes in all of the following discussion. We then generated a grid of halo star models using the 
same input physics; we compare the rotation velocities of Pleiades stars to those for halo composition stars 
with the same disk lifetimes in Figure lb. The distribution of initial conditions will be important for the 
remainder of our results, so we begin by discussing the properties of the Pleiades sample first. We then 
compare the angular momentum properties of our halo and solar abundance models. 



3.1.1. Distribution of Initial Conditions 

The measured rotation rates in the Pleiades are subject to ambiguity because we measure wsini 
rather than the true rotation rate; the true rotation velocities will therefore be higher on average than the 
wsini values. In addition, there are a significant number of upper limits in the data set; this technique 
would therefore miss any population of very slow rotators. Such stars, if present, would have less lithium 
depletion than our technique would permit, and they would therefore influence the predicted range in 
lithium depletion factors. Finally, there is the possibility that the distribution of initial conditions could 
depend on environment; i.e., the assumption that halo stars had the same distribution of disk lifetimes as 
the Pleiades, or that the Pleiades is the same as other open clusters, may be incorrect. 

A distribution of v sin i values can in principle be inverted to infer the true distribution of velocities 



( Chandrasekhar & Munch 195C), but the technique requires a large sample and we would lose the 
information on the mass dependence of rotation which is important for our purposes. Although sini can be 
small, it will be close to 1 for the large majority of stars; only 1/7 of the sample will have sini < 0.5 for a 
random distribution of inclination angles. Furthermore, the young open clusters have a small population of 
rapid rotators and a large population of slow rotators; 13 out of the 91 stars in the Pleiades sample have 
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vsmi > 40 km/s, and 7 have vsini between 20 and 40 km/s. Therefore, only a small fraction of the low 
v sin i systems are intrinsically rapid rotators seen pole-on. For these reasons we have therefore chosen to 
apply an average statistical correction to the data, and do not believe that the shape of the distribution in 
Figure 1 will be greatly altered by inclination angle effects. Rotation periods are not affected by sini, and 
they provide an alternate means of estimating the distribution of initial conditions. Although the present 
rotation period samples are much smaller than the vshii samples, the number of stars with measured 
rotation periods is growing and it may soon be possible to use them in well-studied systems such as the 
Pleiades. 

The issue of the lower envelope of rotation is a more serious concern; there is a bias in the v sin i sample 
against detecting truly slow rotators. We have a series of reasons for believing that the absence of very slow 
rotators is real. First, stellar rotation periods can be measured for stars without the uncertainty introduced 
by the inclination angle; much slower rotation rates can be measured with this technique than with vsini 



data. Long rotation periods are not found in the Pleiades ( Krishnamurthi et al. 1997t) and are very rare 



in the T Tauri population (Edwards et al. 1993, Bouvier et al. 1993), which suggests that any population 



of very slow rotators is small. Furthermore, even disks that last to the main sequence would only produce 
stars with rotation periods of order 10 days, very close to the rotation period corresponding to the vsini 
limits. We would therefore require stars which are locked to their disks with a period much longer than 
that seen in T Tauri stars and with a long disk lifetime to produce a population of ultra-slow rotators large 
enough to affect our results. Finally, progressively younger open clusters have a smaller number of upper 
limits than older systems (there are fewer in a Per than in the Pleiades), and the lower envelope of rotation 
is resolved in the youngest systems (IC 2391 and IC 2602). Recent studies of rotation in the Pleiades are 



also consistent with this conclusion (Bouvier, Forestini, fc Allain 1997) 



The universality of the distribution in Figure la is another important question. KPBS found that the 
distributions of disk lifetimes needed for the different young open clusters were consistent with one another, 
in the sense that the slow rotator distribution peaked at similar disk lifetimes for systems of different ages. 
This does not guarantee that the same will hold true for systems that formed under different conditions, 
such as halo stars and globular cluster stars. The distribution of initial rotation rates cannot be directly 
inferred in old systems because the surface rotation rates of stars converge as they age and lose angular 
momentum. Even in a system as young as the Hyades (600 Myr) the range in surface rotation rates is too 
small to permit a reliable estimate of the initial angular momentum of any given star. It is, however, at 
least a plausible starting point to ask whether the distribution of lithium abundances in different samples 
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is consistent with the distribution of rotation rates seen in young open clusters such as the Pleiades. One 
can in principle invert the question and ask what distribution of disk lifetimes is needed to produce a given 
distribution of lithium abundances. Lithium depletion depends on the rotational history of a given star, 
and therefore the distribution of surface lithium abundances in old stars can be used to obtain information 
about the distribution of initial conditions in systems that are too old to provide direct information on 
the rotation rates. We will show that both the halo star and open cluster lithium depletion patterns are 
consistent with lithium depletion from rotational mixing and with being drawn from the same set of initial 
conditions. 

3.1.2. Metallicity Effects 

For each Pleiades star we can infer a disk lifetime; in Figure lb we compare the rotation properties 
of halo star models with the same age, effective temperature, disk lifetime, and rotation physics (angular 
momentum loss law and angular momentum tranport prescription). We did not run halo star models with 
effective temperatures as low as the bulk of the Pleiades sample, since these stars would not be on the halo 
lithium plateau; we have therefore only shown the T e g range where the two samples overlap. There are 
some important features of this comparison which will have implications for the lithium depletion pattern 
in open cluster and halo stars. First, the overall level of rotation is lower in the halo stars; second, the 
range in rotation rates is smaller. Both phenomena can be traced to differences in stellar structure, and the 
relative pattern is the same as described in PDD and in Chaboyer & Demarque (1994). 

At a given T c ff a halo star will have a lower total mass and moment of inertia; for the same amount 
of angular momentum loss the halo star will therefore spin down more rapidly. Rapid rotators lose more 
angular momentum than slow rotators, which causes their surface rotation rates to converge at late ages. In 
addition, the convergence to similar surface rotation rates occurs more quickly in the halo stars because of 
their smaller moments of inertia. The diffusion coefficients for mixing depend on both the absolute rotation 
rates and the internal angular velocity gradients; both are smaller in halo stars than in solar composition 
stars of the same age and T c ff. The overall level of lithium depletion will therefore be lower in the halo 
star models. In addition, a range in surface rotation rates will produce a range in the degree of rotational 
mixing, and therefore a dispersion in surface lithium abundances. The smaller range in surface rotation 
rates for the halo star models will therefore produce a smaller dispersion in surface lithium abundance at 
fixed T ff. Similar initial conditions will therefore produce different apparent lithium depletion patterns; 
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this is important because such differences are cieariy seen in the observational data. 

3.2. Monte Carlo Simulations of Open Cluster Lithium Depletion 

Our theoretical models can be used to predict lithium depletion in open cluster and halo stars. As 
discussed above, the observed distribution of rotation rates in young open clusters (like the Pleiades) is 
mapped onto a distribution of initial rotation rates by use of an accretion disk model which, during the T 
Tauri phase of its pre-main sequence evolution, locks the rotation period of the star at 10 days. A range in 
decoupling times for these accretion disks will then generate a range of rotation rates for stars on the main 
sequence as a function of time. The rotational properties of the models are calibrated on the rotational 
properties of the open cluster stars and solar models are required to reproduce the solar rotation rate at the 
age of the Sun. This step is independent of the initial conditions for the Sun because different disk lifetimes 
for solar models converge to identical surface rotation rates at the age of the Sun. 

We then need a means of relating the diffusion coefficients for mixing to those for angular momentum; 
the usual approach is to require that a solar model reproduce the solar surface lithium abundance at the age 
of the Sun. Because models with different disk lifetimes have different rotation histories, this step depends 
critically on the disk lifetime assigned to the Sun. If the Sun is assumed to have a typical rotation history 
then solar analogs will have lithium depletion comparable to that in the Sun. On the other hand, if the Sun 
was an unusually rapid rotator (i.e., had a shorter disk lifetime than most stars) the majority of other stars 
will have less lithium depletion. We therefore ran solar calibrated models with the Sun having disk lifetimes 
of (anomalous), 0.3, and 1 ("typical") Myr; these cases will be referred to as sO, s0.3, and si for the 
remainder of the text. The open cluster lithium data are fit poorly when the Sun is assumed to have a disk 
lifetime longer than 1 Myr (see below) . All of our solar disk lifetimes are shorter than the 3 Myr required for 
typical stars; this implies that the Sun was a faster than average rotator in its youth and thus its rotation 
history was not "typical" , strictly speaking. The solar-calibrated models can then be used to predict the 
lithium depletion factor as a function of time for a given disk lifetime, mass, and composition. We then 
investigated the distribution of lithium depletion factors that would be produced by a given distribution of 
disk lifetimes. For each star we chose an initial condition at random from the Pleiades distribution of disk 
lifetimes; we show sample Monte Carlo runs for the four open clusters using two different solar calibrations 
(sO and si) in Figures 2 and 3 respectively. For each cluster all stars were assumed to start with the same 
initial lithium abundance: we chose [Li] = 3.4 for the youngest open clusters and [Li] = 3.3 for the Sun and 
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for M67. Standard model predictions for Li in these clusters are shown as solid curves in Figures 2 and 
3. Note that some main sequence lithium depletion is demanded by the data and that the overall lithium 
depletion pattern as a function of age is in reasonably good agreement with the Monte Carlo simulations of 
the models with rotational mixing. Variations in the initial rotation rates provide a natural explanation for 
dispersion in stellar surface Li/H at fixed T e ff. 

The predictions of the time dependence of lithium depletion due to rotational mixing matches the 
open cluster data best for the case where the Sun was an unusually rapid rotator (Figure 2). This sO 
case corresponds to minimal lithium depletion and predicts modest dispersion in the open cluster lithium 
abundances. As the Sun is made more typical (i.e., as its disk lifetime increases and ZAMS rotation rate 
drops) the overall depletion and dispersion become larger and the overall agreement with the data worsens. 
The intermediate-aged systems are the most sensitive discriminant; the model predictions are similar in 
young systems and the old clusters have less data (which is also less accurate). A more detailed analysis 



of the open cluster lithium depletion pattern will be presented in a paper in preparation ( Narayanan ei 



al. 1998); for the purposes of this paper, we note that open cluster models constructed using the same 
technique that we have applied to halo stars are in overall agreement with the observed lithium depletion 
pattern. 

3.3. Halo Stars 

3.3.1. Parameter Variations 

We can extend the PNK models, designed for rotational mixing in open clusters, to study lithium 
depletion in halo stars. Unlike the open clusters, the halo stars have a range in metal abundance and age. 
We therefore begin with a discussion of the dependence of the lithium depletion pattern on the important 
components of our models: disk lifetimes, solar calibration, metal abundance, and age. In Figure 4 we show 
the effects of changes in all of those ingredients. In all panels the solid line is a base case which corresponds 
to models with a 3 Myr disk lifetime, [Fe/H] = —2.3, and an age of 12 Gyr. 

1. Disk lifetimes. 

A range in disk lifetime can produce significant variations in the predicted lithium depletion, especially 
for the rapid rotators. This will manifest itself as a dispersion in lithium abundance at fixed T c g (top panel 
of Figure 4). However, the large majority of stars in the Pleiades are slow rotators with disk lifetimes 
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greater than 1 Myr. As we will see below, the intrinsic range in depletion factors for slow rotators is very 
small; i.e., only the top 20% of the distribution will have depletion factors greater than or equal to the 
depletion factor of the 1 Myr case see Fig. la), while the remainder will have depletion factors similar to 
that for the 3 Myr case. The overall pattern that would be expected from the Pleiades distribution of initial 
conditions is therefore a majority of stars with a small intrinsic range in abundance with a subpopulation 
(corresponding to rapid rotators) which is overdepleted with respect to the bulk of the stars. 

2. Solar Calibration. 

The solar calibration sets the overall level of depletion in the models, and it will therefore have a strong 
influence on the derived primordial lithium abundance (the second panel of Figure 4). Once the solar 
calibration is set by the open cluster data, it should be the same for all stars ( |Chaboyer and Zahn 1992 ), 
and in particular for halo stars. That is, we expect that the best fit to the halo star lithium abundances 
will correspond to the limit in which the Sun is over-depleted in lithium relative to its counterparts (the sO 
case). 

We define the depletion factor D7 as the ratio of the initial ( 7 Li/H) to the current surface ( 7 Li/H). 
The halo star models of PDD predicted a value of D? close to 10 for plateau stars, with little dependence 
on the model parameters. This behavior is very different from that illustrated in Figure 4, and this can be 
traced to differences in the initial conditions and angular momentum loss law. PDD began with a modest 
range of initial angular momenta (a factor of 10) and applied a loss law which caused a rapid convergence 
in the surface rotation rates. We begin with a wider range of initial angular momenta (disk lifetimes from 
Myr to 10 Myr will produce a factor of 16 range in angular momentum at an age of 10 Myr). The loss 
law that is required by the rotation data also permits a wide range of rotation rates to persist for a longer 
period of time than it did in the models of PDD. The net result is a much larger possible range of lithium 
depletion factors. The models of PDD, for example, had a total range in log(Dy) of 0.25 dex at 0.7Mq and 
[Fe/H]=— 2.3, while the corresponding models in the current study have a range in depletion factors of 0.50 
dex (sO) to 1.30 dex (si). Because the possible range of lithium depletion factors is greatly increased, the 
choice of the initial conditions for the Sun becomes correspondingly more important. In effect, the earlier 
generation of models had a small enough effective range in initial conditions that the Sun was by definition 
a typical star and the typical depletion factors were therefore forced to be similar to those for the Sun. If 
the Sun is treated as a typical star with a disk lifetime of 3 Myr, we recover the mean depletion factor of 
PDD (log(Z?7) = —1.0) with a much larger range of depletion factors than permitted by the data. 
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3. Metal Abundance 

Models with different metal abundance differ in their overall lithium depletion factors by a modest 
amount, less than 0.1 dex change in lithium depletion per 1.0 dex change in [Fe/H]. The overall trend is in 
the sense that the more metal-rich stars deplete more than the more metal-poor stars; this is to be expected 
in light of the relative differences between solar composition and halo star rotational properties illustrated 
in Figure 1. 

This is opposite to the trend found by Thorburn (1994) in the observational data, for which stars with 
lower metal abundance have, on average, slightly lower lithium abundances, by 0.1 dex in [Li] for 1.0 dex 
in [Fe/H]. However, there are some possible complications in this relative comparison. The slope of the 
[Li]- [Fe/H] relationship depends upon the assumed model atmospheres. Furthermore, the derived value of 
the slope depends on the data analysis method which is used and is different for different subsets of the 



data ( Ryan et al. 1996 ). By contrast, the [Li]-T e ff slope reported by Ryan et al. (1996) is resistant to these 
effects. 

Within the framework of the rotational mixing hypothesis, one would also expect stochastic variations 
in the disk lifetimes and therefore an additional random component in the lithium abundances. Since about 
20% of the Pleiades stars have modest to rapid rotation, the sample could (for example) happen to have 
slighty more rapid rotators at low [Fe/H] and slightly fewer at relatively high [Fe/H] by chance; this effect 
is not included in the Ryan et al. (1996) or Thorburn (1994) analysis. 

Another possibility is that the trends in the models and the data are correct, and that Galactic 
production of 7 Li is responsible for the difference between the trends. In all of the models we require some 
lithium production during the lifetime of the Galaxy; this would imply that the more metal-rich stars 
begin with slightly more lithium, which would counteract the [Li]-T c g trend predicted by the models. The 
Galactic evolution of lithium remains a mystery, but we can estimate the effect as follows. If we scale 
the solar system lithium abundance ([Li] =3.3) linearly with metallicity the newly produced lithium will 
contribute A [Li] = 0.0, 1.0, 2.0 at [Fe/H] = -3.3, -2.3, -1.3 respectively Note that if the post-BBN 
lithium production scales faster than linearly with metallicity such halo contamination will be negligible. 
We infer a primordial abundance in the range [Li] = 2.35 — 2.75. The total abundance at the low end of the 
range (primordial plus production) would therefore be 2.35, 2.37, and 2.51 at [Fe/H] = —3.3,-2.3,-1.3 
respectively; by comparison, at the high end of the range the comparable abundances would be 2.75, 
2.76, and 2.82. This suggests that galactic production could influence the Li- [Fe/H] slope at higher metal 
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abundances but not for abundances below —2.3. 

As we will show below, the detection of 6 Li in HD 84937 requires some 6 Li and Li production at early 
epochs; this production could have been highly variable and possibly even largely absent for the most metal 
poor stars. We therefore cannot exclude an additional production component between the very metal-poor 
stars and the more typical plateau stars, but caution that there may not be a reliable relationship between 
production and metal abundance in this case. If the level of 6 Li seen in HD 84937 is typical for stars at 
[Fe/H]=— 2.3, and the production of 7 Li by a — a fusion is comparable to that of 6 Li, then we can set an 
upper limit on the variation in the initial 7 Li abundance between the most metal poor stars and the typical 
plateau stars by computing the 7 Li yield required to produce the observed 6 Li in HD 84937. The 7 Li yield 
required to explain the observed 6 Li in HD 84937 is in the range of 3 xlCP 11 to 7 xlO~ n . This suggest 
that variable production of 7 Li could have a modest impact on the [Li]-[Fe/H] slope even at early epochs. 

In light of the uncertainties discussed above we have adopted an empirical approach. The value of 
the [Li]-[Fe/H] slope which minimizes the dispersion in the data set is computed, and the observed lithium 
abundances are then corrected for the effects of mctallicity. We have therefore used the [Fe/H] — —2.3 
models for our theoretical calculations and corrected the Thorburn (1994) abundances to a uniform [Fe/H] 
at this level using the relationship of Thorburn (1994), 

Mdctrcndcd = [^observed + 0.l3([Fe/H] + 2.3). (4) 

This raw data is compared with the data corrected to a uniform abundance of —2.3 in Figure 5. We 
note that another approach would be to draw the distribution of theoretical depletion factors from the 
distribution of [Fc/H] in the Thorburn sample, i.e., to Monte Carlo both the initial conditions and the 
[Fe/H] abundances. In numerical tests the correction in equation (4) reduces the dispersion in the data 
set by an amount comparable to the level introduced by the observed distribution of metal abundances 
convolved with the theoretical models. 

4. Age. 

The final ingredient is age; as shown in the bottom panel of Figure 4 even significant variations in age 
produce only small changes in the lithium abundance at fixed T e g. This results from a combination of a 
low rate of lithium depletion at old ages and the evolutionary changes in T e g as a function of time, which 
move models along the plateau. We therefore neglect the impact of age variations on the dispersion and 
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depletion patterns. 

3.3.2. Monte Carlo Simulations: Trends with T e ff md Dispersion 

Monte Carlo simulations of the lithium abundances of halo stars are compared with the Thorburn 
(1994) dataset detrended to a uniform [Fe/H] of —2.3 in Figure 6. The different models all show a [Li]-T e ff 
trend similar to the data. Previous models with rotational mixing had a downward curvature at the hot end 
of the plateau, which is highly suppressed in these models, although there is a slight downward curvature 
in the si case. This difference can be attributed to improvements in the angular momentum loss law. The 
open cluster rotation data requires stronger angular momentum loss for higher mass stars; this acts to 
preferentially suppress lithium depletion in the higher mass models at the hot end of the plateau relative 
to models for cooler lower mass stars. 

The Pleiades distribution of initial conditions also produces a mean trend with small internal scatter 
and a smaller population of overdepleted stars. As the solar disk lifetime is increased the dispersion and 
overall depletion both increase; we can therefore use the dispersion in the halo plateau to set constraints 
on the overall depletion. The simulations in this section do not include observational errors, whose impact 
on dispersion estimates will be discussed in the next section. The Thorburn sample also contains a small 
population of highly overdepleted stars; their existence is a problem for standard stellar models and they are 
usually ignored in discussions of the dispersion in the halo plateau. We will not ignore the outliers; in fact, 
such objects are a natural consequence of the models with rotational mixing combined with the observed 
distribution of rotation rates in young clusters. The maximum depletion level for the rapid rotators also 
increases as the absolute depletion increases. The existence of highly overdepleted stars can therefore be 
used to constrain the minimum level of absolute depletion in the halo plateau; if the absolute depletion in 
the halo plateau is too small then stars far below the plateau should not be present, but they are. 

Both the trends in the [Li]-T e ff plane and the dispersion in the plateau have been used to argue 
against significant lithium depletion in halo stars. The simulations in Figure 6 indicate that the current 
class of models are compatible with both, but the agreement between the data and the models worsens 
if the absolute depletion is too high. Similar trends are present in the open cluster simulations. In other 
words, the open cluster and halo star data are consistent with being drawn from a ZAMS angular momenta 
distribution similar to that required to explain the distribution of rotation rates in young open clusters. 
Furthermore, both are best explained by a solar calibration in which the Sun had an accretion disk lifetime 
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shorter than the typical value for low mass stars. 

In the next section we quantify the constraints placed on the overall lithium depletion by the dispersion 
and the highly overdepleted stars. We also discuss the relative lithium-6 and lithium- 7 depletion factors 
and constraints on the overall destruction of lithium-7 from the detection of lithium-6 in HD 84937. 

4. Constraints on Rotational Mixing from Halo Stars 
4.1. Dispersion of the Plateau Abundances 

Lithium abundances versus effective temperatures for warm (T e g > 5800K), metal-poor ([Fe/H] 



< —1.3) stars (as observed by Thorburn 1994 ) were shown in Figure 5. Four independent groups have 



observed Popn plateau stars (Spite fe Spite 1982; Spite et al. 1984; Rebolo et al. 1988; Thorburn 1994) 



The differences in the lithium abundances derived by different observers from Popn halo stars can be 
traced to the choice of temperature scale and of model atmosphere required to extract an abundance from 
an observed absorption line. A survey of the literature shows that there may be a systematic uncertainty 
in the lithium abundance scale of order ±0.1 dex due to the choice of atmosphere model and of the stellar 
temperature scale. The intrinsic star-to-star dispersion is likely to be much smaller than this. There may 
also be correlations in the data in the sense that the lithium abundance increases slightly with effective 
temperature and with metallicity (Thorburn 1994; Ryan et al. 1996] ) or not (Molaro, Primas, & Bonifacio 



1995| ). We adopt [Li]p op n = 2.25 ± 0.10 as the weighted mean plus estimated systematic error of all the 



plateau data. The raw dispersion about the trend is 0.16 dex (Thorburn 1994). When the T e g and [Fe/H] 
trends are removed, Thorburn estimated the dispersion as 0.13 dex; by comparison, her estimate of the 
observational error is 0.09 dex. Estimates of the intrinsic star-to-star dispersion are therefore roughly half 



the raw dispersion ( 0.08 dex (Deliyannis, Pinsonneault, fe Duncan 1993), 0.06 dex (Thorburn 1994), < 0.1 



dex (|Bonifacio fc Molaro 1997[ )) 



In our rotationally mixed models, the "anomalous Sun" case corresponds to lithium depletion of 0.2 
dex and a dispersion in depletion of 0.09 dex, consistent with estimates of the intrinsic dispersion. The 
intermediate disk case corresponds to lithium depletion of 0.4 dex and a dispersion in depletion of 0.16 
dex which is at the limit of the raw dispersion. The "typical Sun" case corresponds to lithium depletion 
of 0.69 dex and a dispersion of 0.25 dex which is considerably larger than the raw dispersion. This initial 
comparison suggests that an upper-bound to lithium depletion in halo stars is 0.4 dex. If the Sun is 
anomalous (as favored by the open cluster data) then little dispersion in the plateau is predicted, consistent 
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with the data, along with a smaller overall depletion of 0.2 dex. So, two indicators, the open cluster data 
along with the lack of significant dispersion in the halo lithium abundances, both point to halo star lithium 
depletion in the range 0.2 - 0.4 dex. 

However, the calculation of the observational dispersion excludes 6 upper limits, and the distribution of 
theoretical depletion factors is significantly non-gaussian. We can use these additional characteristics of the 
plateau data to constrain the absolute depletion. Standard models would require that the only dispersion 
in halo star lithium abundances is due to age, metal abundance, and observational error. Age is a small 
effect in all classes of models, and the effects of metal abundance can be included by detrending the data. 
We would therefore expect the shape of the distribution to be determined by the observational errors and 
thus it should be gaussian. 

By contrast, the distribution of lithium depletion factors in rotational models is skewed. The majority 
of young cluster stars are slow rotators, so the majority of the depletion factors are tightly clustered with a 
very small intrinsic dispersion; there is a component of stars with moderate rotation rates which produce a 
population with larger depletion factors and a small fraction (2-3%) of rapid rotators with large depletion 
factors. In addition, the distribution will be smeared out by observational errors. We ran 1000 Monte Carlo 
simulations at intervals of 50K across the plateau and drew theoretical lithium depiction factors from the 
Pleiades distribution of initial conditions. In one set of runs observational errors were neglected; in a second 
set of runs we added a random error of 0.09 dex. The lithium depiction factors for each T c ff were then 
sorted. In Figure 7 we compare the detrended Thorburn data (bottom panel) with the 2.5%, median, and 
97.5% bounds of the theoretical depletion factors for our sO case (top), s0.3 case (second panel) and si case 
(third panel). The zero point was set by requiring that half the stars in the sample for all three cases were 
above the median and half the stars were below it. Figure 8 presents the same cases for simulations which 
include observational errors. The top and bottom dashed lines can be regarded as effective ±2cr bounds. 
Note that the observed lithium abundances have a similar morphology to the simulations; in particular, 
they are not gaussian distributed about the mean as would be expected if these stars had suffered no lithium 
depletion. Furthermore, the intrinsic dispersion for the majority of stars is significantly less than the 
observational errors; in other words, the dispersion for the majority of stars is dominated by observational 
errors. The theoretical models predict that there should be a small number of highly ovcrdcpleted stars, 
which are indeed observed. 

This can be further quantified by comparing the distribution of abundances about the median for the 
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simulations to the observed distribution of abundances about the median. We chose a reference temperature 
of 6000K and an [Fe/H] of —2.3 for the simulations and used the Thorburn (1994) mean trend to correct 
the data to a uniform T e ff and [Fe/H]. Histogram plots of the number of stars as a function of lithium 
abundance are compared to the number of stars as a function of lithium depletion factor expected in 
the simulations without and with observational errors in Figures 9 and 10 respectively. The sO case is a 
good match to the width of the distribution, while the s0.3 case is in marginal agreement and the si case 
clearly has a dispersion too large to be compatible with the data. For comparison purposes, an analogous 
distribution of depletion factors for solar abundance models at a temperature of 6000K is also shown. 

The dispersion in the bulk of the halo population is therefore consistent with overall depletions in 
the range of 0.2-0.4 dex, with lower depletion factors being favored. However, the existence of highly 
overdepleted stars can also be used to set bounds on the absolute depletion, and this favors higher absolute 
depletion factors in the above range. For the sO case, the median depletion is 0.2 dex and the maximum is 
0.6 dex; stars more than 0.4 dex below the plateau should therefore not be seen, in contradiction with the 
data. The corresponding bounds for the s0.3 case are 0.4 and 1.2 dex, permitting stars as much as 0.8 dex 
below the plateau; this is in better agreement with the data. We also note that the observed fraction of 
such stars (4 out of 90) is similar to the fraction of very rapid rotators in young open clusters. We therefore 
conclude that the distribution of lithium abundances in the plateau is consistent with a range of 0.2-0.4 
dex absolute depletion, with the bulk dispersion favoring lower values and the overdepleted stars favoring 
higher values within this range. 



4.2. 6 Li Depletion 

In Figure 11 we show the relationship between 7 Li destruction and 6 Li destruction in the rotationally 
mixed models chosen to correspond to the hot halo subgiant HD 84937 ([Fe/H]« —2.2 and 6 Li/ 7 Li w 0.06) 
( |Smith et al. 1993| |Hobbs k Thorburn 1994|jl997| ), to date the only star in which the detection of 6 Li 



has been confirmed by more than one group (see Hobbs & Thorburn (1997) for a detailed observational 
summary of the 6 Li data). Note that, depending on the solar normalization, the ratio of 6 Li destruction to 
7 Li destruction can vary from a few (for the sO case) to 10 (for the si case). The absolute 6 Li depletion 
factor varies between roughly 6 and 70, respectively. Note that this is not the usual folklore of lithium 
depletion which anticipates substantial 6 Li destruction in any star where even minimal 7 Li destruction has 



occurred (Brown & Schramm 198S , Vauclair & Charbonnel (1998)). Indeed, it has often been argued that 
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because 6 Li has been observed in at least one halo star, there can be little, if any, 7 Li depletion in any 



Pop ii plateau stars (Steigman et al. 1993). However, this argument was based on standard models without 
rotational mixing where all the lithium destruction is by convective burning during the star's approach to 
the main sequence. Then, since 6 Li is much more fragile than 7 Li, any star in which 6 Li is not significantly 
destroyed is not expected to destroy any 7 Li at all. This 6 Li - 7 Li connection is relaxed in models which 
include rotational mixing. In the convective burning only models, all the material on the surface has been 
exposed to the same temperature whereas in models with rotational mixing different fractions of the surface 
material are exposed to different temperatures. Although the rank ordering (with binding energy) of the 
magnitude of lithium destruction is preserved (i.e., 6 Li is always destroyed more than 7 Li), the relative 
amounts of depletion can vary since the surface of the star is a mixture of gas which has been exposed 
to different temperatures. Note also that rotational models predict a dispersion of 6 Li depletion at fixed 
T e ff corresponding to different initial rotational velocities as drawn from a distribution of initial conditions. 
Below we use the Popn 6 Li observations to further constrain models of rotational mixing and overall 
lithium depletion. 

To use the observed abundance of 6 Li as a constraint on the amount of lithium depletion, we require 
an estimate of its initial abundance. 6 Li is only produced in interesting (i.e., observable) amounts by 
cosmic ray nucleosynthesis via two channels: (1) cosmic ray proton and alpha spallation of CNO nuclei 
in the ambient ISM (and vice-versa) and, (2) alpha-alpha fusion reactions. To avoid introducing more 
adjustable parameters than can be constrained by the observational data we can tie the CNO contribution 
of 6 Li to the observed abundances of B and Be which are produced by the same spallation process. The 



Be and B abundances in Popn halo stars are observed to increase linearly with Fe/H (Duncan et al. 1997) 
and, therefore, so too should the "CNO-tagged" 6 Li. The predicted ratio of spallation-produced 6 Li to 
Be or B depends on the specific cosmic ray spectrum and relative target/projectile abundances adopted. 
Some of these details can be avoided if we simply normalize 6 Li production to the solar abundance of 6 Li 
([ 6 Li] Q ~ 2.2): [ 6 Li] CR(CNO) = [ 6 Li] + [Fe/H] = 2.2 + [Fe/H]. This represents an upper limit to the CNO 
CR production of 6 Li at any epoch. Since the metallicity of the Spite plateau stars is [Fe/H] < —1.3, the 
CNO contribution will, in general, be quite small: [ 6 Li]cr(cno) < 0.9. If this were the only production 
mechanism for 6 Li in halo stars, it would, along with the plateau value of [Li] ~ 2.2, correspond to a ratio 
6 Li/ 7 Li < 1/20, which is comparable to that observed. In this case (production via spallation of CNO) 
there would be room for little, if any, 6 Li destruction. 

To quantify this point, consider the case of HD 84937, a hot (T c g ~ 6300K) metal-poor star in which 
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Li observations have been claimed by two independent groups ( Smith et al. 1993| ; Hobbs & Thorburn 
1994], 1997). They find [Fc/H] = -2.2 ± 0.2 and [ 6 Li] OB s = l-Otoi, an order 01 magnitude more 6 Li than 



the maximum predicted from CNO spallation alone: [ 6 Li]cpj(cNo) = 0.0 ± 0.2. It must be emphasized 
that this can NOT be used as evidence for the absence of 6 Li destruction in halo stars but instead points 
to the necessity of an additional source for the observed 6 Li (in particular, one which does not scale with 
metallicity) . 



Alpha- alpha fusion synthesis of 6 Li ( Steigman fc Walker 1992 ; SW), is not necessarily coupled to 
the cosmic ray spallation production of Be and B and could, in principle, be such a source of 6 Li (and a 
comparable amount of Li). Naively we expect the fusion contribution to scale as one lower power of Fe/H 
than does the CNO-spallation contribution (i.e., that it may be nearly independent of metallicity). As 
before for spallation, the normalization of the alpha-alpha 6 Li contribution is essentially arbitrary. As an 
extreme upper limit to the abundance of alpha-alpha produced 6 Li in the gas out of which the Popn stars 
form, we may adopt the solar abundance: 



l 6 Li} C R {aa ) < 2.2. (5) 

This bound is certainly robust (and most likely overly conservative) since it assumes that CR fusion in 
the very early Galaxy generated 6 Li at the solar abundance. Nevertheless, it represents the maximum 
amount of 6 Li any Pop II star could contain. Defining the 6 Li depletion factor, D§, as: ( 6 Li/H)oss = 
{ 6 U/U) CR{aa) /D 6 , we have for HD 84937 

logA, < 1.2±g;f. (6) 

That is, there is a model indenpendent upper bound to 6 Li destruction of a factor of 14-25 for HD 84937. 
The rotational models for HD 84937 which violate this constraint (i.e., destroy too much 6 Li) are the same 
ones which generate too much dispersion in the plateau abundance of 7 Li and do not correctly predict the 
behavior of the lithium abundances in open clusters. This Dq constraint allows us to obtain an independent 
upper bound on the depletion of 7 Li in Popn plateau stars (see Figure 11) of 0.5 — 0.6 dex. We comment 
here that although Lemoine et al. (1997) claim to have constrained D§ to be strictly less than four, they 
did not include the fusion production of 6 Li via alpha-alpha which we have seen may be required if 6 Li has 
been detected in HD 84937. 
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4.3. Relative Depletion of 6 Li and 7 Li 

Since post-BBN production of lithium can compete with surface depletion, lithium production and 
destruction are coupled. Thus the dispersion in the Popn lithium data itself can be used to bound any 
contributions to the observed Popn lithium abundance from post-big bang production and/or stellar 
destruction. Every Popn star should start with the same amount of BBN-produccd 7 Li and variations 
in post-BBN production and stellar destruction will contribute to dispersion in the observed "plateau" 
abundances. Any observational limits to the dispersion provide a limit to the combination of CR-production 
and stellar depletion of lithium in the halo stars. That is, for any halo star, the observed total ( 7 Li plus 
6 Li) lithium abundance can be expressed as 

OBS _ 2/7BBN + J/7CR 2/6CR 

Vlz — p; 1 pr - ! {') 

where the various j/7's represent the observed (OBS), BBN, and cosmic ray produced (CR) number fractions 
of 7 Li relative to hydrogen. The CR contribution to the observed lithium, with account taken of stellar 
depletion, may be bounded by the dispersion in the observed lithium abundances: 



2/7CR , V6CR 



< Ay2f S « 0.8 x 1(T 10 , (8) 



D 7 D 6 

where we have assumed that the BBN contribution is fixed by the minimum of the plateau: 

^= = (y OBS )mm „ lA x 1Q -10 (Q) 

To bound D 6 /D 7 we may then write for the (inverse of the) 6 Li fraction 



V Li )oBS D 7 \vJcR 



j/7BBN ^ 
V7CR 



1. (10) 



These relations can be used along with the observational data to bound the combinations 2/7BBN/-D7 and 
2/7CR/-D7, respectively. The 6-to-7 ratio expected (SW) from CR nucleosynthesis is (ye/V7)cR < 0.9, leading 
to a bound on Dq/Dj: 

^<8. (11) 

That is, the small dispersion of the Pop 11 lithium data along with the 6 Li abundance limits the relative 
6-to-7 depletion. We caution that this calculation relies on the measured 6 Li in only one star, HD 84937. 
We find that the same models with rotation which, earlier, failed to reproduce the open cluster lithium 
abundance patterns, that predicted too much dispersion in the Pop 11 lithium plateau, and too much 
absolute destruction of 6 Li, also violate this constraint on the relative depletions of 6 Li and 7 Li. From 
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all these - independent - constraints a consistent picture of non-zero but limited lithium depletion in the 
warm, metal-poor Popn plateau stars emerges. 

In conclusion, we find that those models with rotational mixing which are constrained to fit the lithium 
abundance patterns in open clusters require that the Sun be an anomalously fast early rotator which has 
overdepleted lithium relative to most stars of its mass and age. Extending these same models to halo stars 
yields depletion patterns which can be constrained by the lithium abundances observed in these stars. The 
predicted dispersion in lithium abundances and the predicted 6 Li depletion factors are consistent with the 
Pop ii data for the same solar normalization required by the open cluster data. Thus, models which satisfy 
all of the above criteria suggest that halo stars must deplete their initial lithium by at least 0.2 dex, but 
by no more than 0.4 dex. The conclusion that halo stars must deplete lithium in this range relies only on 
the assumption that halo stars have the same distribution of initial rotational velocities as do the open 
cluster stars. If halo stars do not rotate then standard models would predict no lithium depletion except 
at low T e ff, contrary to the data. Alternatively, significantly more lithium might be destroyed without 
increasing the dispersion in the observed lithium abundances if every halo star had nearly identical angular 
momentum histories. 

Armed with our estimate of the lithium depletion factor relevant for the Spite plateau halo stars we 
may use the data to bound the primordial abundance of lithium. 



Adopting for the Spite plateau abundance [Li]oBS = 2.25 ± 0.10 and for the depletion factor log 
D 7 = 0.30 ± 0.10, the primordial abundance should lie within the full range 2.35 < [Li]p < 2.75 or, 



Here we have assumed the errors in both the observed lithium abundance and the depletion factor are 
uniformly distributed with a width of 0.2 dex. This results in a triangular likelihood distribution for [Li]p 
which has a 95% C.L. range for the log of the primordial lithium abundance only 5% smaller than the full 
range values we quote here. Since the shape of the distributions we assume are, at best, an approximation, 
we quote the full range in [Li] p resulting from the convolution of the two distributions. 



[Li] P = [Li] BS + log£>7 



(12) 



2.2 < 10 w (Li/H) P < 5-6. 



(13) 
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5. Implications for BBN 

The predictions of standard BBN are uniquely determined by one parameter, the density of baryons, 
parameterized by r\ - the baryon-to-photon ratio: 7710 = 273flsh 2 (J1b is the ratio of the baryon density to 
the critical density and the Hubble parameter is Ho = lOOh km/s/Mpc; 7710 = lO 10 ??). For standard BBN 
the early universe is assumed to have been homogeneous and expanding isotropically and the energy density 
at the time nucleosythesis begins (about 1 second after the big bang) is described by the standard model of 
particle physics (ptot = Py + Pe +N l/ p l/ , where p 7 , p e , and p u are the energy density of photons, electrons and 
positrons, and massless neutrinos (one species), respectively, and N v is the equivalent number of massless 
neutrino species which, in standard BBN is exactly 3). In Figure 12 the primordial abundances predicted 
by standard BBN are shown as a function of 77. The width of each curve reflects the ±2tr uncertainty in the 
predictions. 

Our analysis of lithium depletion by rotational mixing, when combined with the Pop 11 lithium data 
leads to a primordial lithium abundance in the range: 2.2 < 10 10 (Li/H)p < 5.6. Comparing with the 
predictions of standard BBN (see Figure 12), lithium identifies two possible ranges for 77. A "low-77" branch: 
0.8 < 7710 < 1.7 and a "high-77" branch: 3.7 < 7710 < 9.0. Once r\ is identified the primordial abundances of 
deuterium and helium-4 (also, helium-3) follow from the predictions of standard BBN. 



5.1. Low-77 



The baryon density corresponding to the low-77 range is quite small: 0.003 < JIb/i 2 < 0.006. If the x-ray 



emission from rich clusters of galaxies provides a "fair" sample of the universal fraction of baryons ( White 
et al. 1993| ; jStei gman fc Felten 1995| ; |Evrard 1997| ; |Stcigman, Hata fc Felten 1997[ ) this argues strongly for a 
very low density, open Universe: £1 = Ob/Zb is O.l/i -1 / 2 . Such a low baryon density may also be in conflict 
with the baryon density required to reproduce the Lyman— a forest data in large-scale structure simulations 
(Weinberg et al. 1997). This low-77 range predicts a high primordial deuterium abundance (10 5 (D/H)p > 
17) which is inconsistent with the low deuterium abundance inferred for some high-redshift, low metallicity 



QSO absorption systems (Buries fe Tytlcr 1996; Tytlcr, Fan, fc Buries 1996; Tytlcr, Buries, fc Kirkman 



1996 ; Buries fc Tytler 1997 ) but not with the high abundances inferred for some other systems ( Songaila et 
al. 1994j |Carswell et al. 1994| ; |Rugers fc Hogan 1996a| ; |Rugers fc Hogan 1996b] ; |Songaila et al. 1997| ; |Webb| 



et al. 1997). This low-77 range corresponds to a BBN-predicted helium abundance (Yp < 0.235) which is in 
excellent agreement with that inferred from emission-line observations of the low-metallicity, extragalactic 
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Hn regions flOlivc fc Stcigman 1996|; plivc, Skillman, fc Stcigman 1997[ |Hogan, Olive, fc Scully 1997]). 



5.2. High-?/ 



Avoiding the minimum in the BBN lithium "valley" (see Figure 12), this alternate branch favors a 
higher baryon density: 0.014 < Q,^h 2 < 0.033 which, when combined with the x-ray cluster baryon fraction 
estimate, permits a higher total density Universe: 0.2 < Qh 1 / 2 < 0.6. A baryon density in this range is 
consistent with the Lyman— a forest data (Weinberg et at 1997). On this branch the primordial deuterium 
abundance is predicted to be smaller (6.6 > 10 5 (D / H)p > 1.1), consistent with some recent estimates from 
high-z QSO absorbers. However, the primordial helium abundance does pose a challenge to standard BBN 
since 0.243 < Yp < 0.253 is predicted while the Hn region data suggest 0.230 < Yp < 0.238 (but see also 



[zotov, Thuan, fc Lipovetsky 1994, 1997 and Thuan and Izotov 1997 who derive 0.244 ± 0.002 from an 



independent data set). 



6. Summary 

The lithium abundances of stars have interesting implications for cosmology and for the theory of 
stellar structure and evolution. The conclusions of this paper impact both of these areas. Standard stellar 
models implicitly assume that a variety of physical processes known to occur in real stars can be neglected; 
for many purposes this is a good approximation. There is now extensive evidence in Popi stars that the 
predictions of standard models are not in agreement with all the data. Furthermore, work from a variety of 
investigators on physical effects not included in the standard models has concluded that there are physically 
well-motivated processes which can affect the surface lithium abundances of stars. Standard models and 
the previous generation of nonstandard models predicted initial 7 Li abundances that differ by up to a factor 
of ten from the same set of current Pop n data, with very different implications for cosmology. 

We believe that mild envelope mixing driven by rotation is the most promising candidate for explaining 
the complete observational picture. An improved treatment of angular momentum evolution in low mass 
stars now allows us to generate stellar models with rotation which are consistent with the rotation rates 
observed as a function of mass and age. We have inferred the distribution of initial angular momenta in 
young clusters and computed the distribution of lithium depletion factors as a function of mass, composition, 
and age. The resulting lithium depletion pattern is in good agreement with both Popi and Popn data, 
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and the distribution of theoretical depletion factors is consistent with the distribution of abundances. 
For Pop ii stars in particular, the observed slope of the [Li]-T e g relationship, the observed dispersion in 
abundance at fixed T e ff, the existence of a small population of overdepleted stars, and the simultaneous 
detection of 6 Li and 7 Li in one halo star are all consistent with the predictions of the theoretical models. 
The primary uncertainty in the absolute depletion of 7 Li is the initial angular momentum of the Sun, which 
is used to calibrate the mixing diffusion coefficients. 

We have used the observed properties mentioned above to set bounds on 7 Li depletion for Popn stars. 
Those models which fit the Popi data best predict a small but real depletion of lithium in the warm 
(T eff > 5800K), metal-poor ([Fe/H] < -1.3) Popn stars: 0.2 < logD 7 < 0.4. These same models are 
consistent with the very small observed dispersion in abundances about the plateau value and with the 
survival of some 6 Li. These depletion factors are significantly less than in previous models with rotational 
mixing (PDD, Chaboyer & Demarque 1994), and this difference can be directly attributed to the adoption 
of an improved treatment of angular momentum evolution. The uncertainties in the 7 Li depletion factor 
can be reduced by a combination of new data and improved stellar evolution models. A large set of 
accurate lithium abundances in old open clusters would enable us to calibrate the diffusion coefficients 
for mixing without relying on the solar calibration. Microscopic diffusion and rotational mixing should be 
considered together, although the work of Chaboyer & Demarque (1994) indicated that the predictions 
of such models are similar to those which include rotation alone. The observed 7 Li depletion pattern is 
relatively insensitive to the treatment of internal angular momentum transport, with a factor of ten in the 
time scale corresponding to less than a 10% change in the logarithmic lithium depletion factor in the work 
of PDD. However, the solar rotation curve does provide evidence for angular momentum transport from 
mechanisms not considered in the current paper, such as magnetic fields and/or internal gravity waves. 
Lithium depletion in models which include these effects and hydrodynamic mechanisms should be explored. 
Intermediate metallicity stars could also provide a valuable test of the dependence of lithium depletion on 
metal abundance. 

When our constraints on lithium-depletion are combined with the abundances inferred from observations 
of the Spite plateau halo stars ([LijoBS = 2.25± 0.10) we find for the range in the abundance of primordial 
lithium: 2.35 < [Li] P < 2.75 (2.2 < 10 10 (Li/H) P < 5.6). Although lithium is not the ideal baryometer, 
a comparison with the predictions of standard BBN identifies two options: low-77 and high-ry. The low-77 
branch suggests a very low baryon density, open Universe which may be in conflict with the baryon density 
inferred from observations of the Lyman-alpha forest. Although the helium abundance predicted for this 
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low-77 option agrees with that inferred from H n region data, the very high predicted deuterium abundance 
may not be consistent with observations. In contrast, the high-?7 branch corresponds to a baryon density 
consistent with the Ly-a data and permits a higher density Universe. In this case the predicted primordial 
deuterium abundance is in excellent agreement with the low deuterium QSO data but a considerably higher 
primordial helium abundance is predicted than is inferred from some of the observational data. 

This work is supported by the Department of Energy Contract No. DE-AC02-76-ER01545 at The Ohio 
State University and NASA ATP Grant 68169492. TPW thanks Sylvie Vauclair, Con Deliyannis, Corinnc 
Charbonnel, Francois Spite, Rafael Rebolo, and Paulo Molaro for useful conversations. 
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Fig. 1. — (a) Distribution of accretion disk lifetimes inferred from the observed rotational velocities of 
Pleiades stars. Data is taken from Soderblom et al. (1993a). (6) Measured rotational velocities of the 
Pleiades stars (crosses) in km/s and the rotational velocities of the halo composition models (solid circles) 
at an age of 70 Myr with the same T c g and disk lifetimes. 

Fig. 2. — Comparison of the [Li] measured in four different open clusters (crosses) to the [Li] predicted by 
the sO models (filled squares) in a typical Monte-Carlo run. The open inverted triangles represent those stars 
for which only an upper limit is quoted in the observations. Lithium data are taken from Soderblom et al. 
(1993b) and Balachandran (1995). The solid lines are standard model 7 Li depletion factors for the metal 
abundances of the clusters ([Fe/H]= —0.03, +0.12, +0.04, and 0.0 for the Pleiades, Hyades, Praesepe, and 
M67 respectively). The initial abundance [Li] in the theoretical models was assumed to be 3.3 for M67 and 
3.4 for the other clusters. 

Fig. 3. — As Figure 2, except for a Monte Carlo simulation using the si case. 

Fig. 4. — Sensitivity of the lithium depletion patterns to changes in (a) the time scale for coupling between 
the accretion disk and the protostar in Myr (r^), (fc) the solar calibration (accretion disk lifetime for the 
Sun), (c) [Fe/H] and (d) age. In all these panels, the solid line represents a base case for a typical halo star 
with Td = 3 Myr, [Fe/H] = —2.3, sO and age = 12 Gyr. D-j is defined as the ratio of the current surface 7 Li 
to the initial 7 Li. 

Fig. 5. — [Li] in halo stars before (panel a) and after (panel b) correcting for the metallicity trends. 
The abundance data were taken from Thorburn (1994) and corrected to a uniform metal abundance using 
equation 4. 

Fig. 6. — Sample Monte-Carlo distributions of 7 Li depletion factors for the three different solar calibrations, 
(a) sO, (b) s0.3 and (c) si. (d) Thorburn (1994) lithium data corrected to a uniform metallicity of [Fe/H] = 
-2.3. 

Fig. 7. — Lithium depletion factors for the three different solar calibrations before accounting for the 
observational errors, (a) sO, (6) s0.3 and (c) si. The solid line shows the median value, while the dashed 
lines show the 2.5 and 97.5 percentile values. The filled circles show the observed [Li] in the halo stars after 
correcting them all to the same metallicity of [Fe/H] = —2.3. 

Fig. 8. — As Figure 7, except a random observational error of 0.09 dex was included in the theoretical 
simulations. 

Fig. 9. — Distribution of the lithium depletion factors for the three different solar calibrations before 
accounting for the observational errors. All the depletion factors are estimated at T c fj=6000K and 
[Fe/H]=— 2.3. The top three panels show the results for the three different solar calibrations, (a) sO, (b) s0.3 
and (c) si. The dashed histograms are the distributions of depletion factors for the same T e g of 6000K but 
with a solar [Fe/H]. Panel (d) shows the observed distribution of [Li] in the halo stars after correcting them 
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all to a metallicity of [Fc/H] = —2.3. A correction of 0.04(T e fj — 6000-fQ was also applied to the data to 
correct for the dependence of [Li] on effective temperature. 

Fig. 10. — As Figure 9, except a random observational error of 0.09 dex was included in the simulations. 

Fig. 11. — 6 Li depletion as a function of 7 Li depletion for models of HD 84937. Three different solar 
calibrations are shown. For each calibration the 6 Li and Li depletion factors are indicated by the symbols 
(squares, open circles, crosses) for disk lifetimes of 1, 2, 3, 5, and 7 Myr (running from lower left to upper 
right respectively). We have connected these discrete points by solid lines in each case. A least squares fit 
to 6 Li depiction as a function of 7 Li depletion for a disk lifetime of 3 Myr is also indicated. 

Fig. 12. — The predicted standard BBN yields (dashed lines represent 2-<r uncertainties with 3 massless 
neutrino species) of light elements ( 4 He mass fraction, D and 7 Li number fraction relative to hydrogen) 
as a function of r], the baryon-to-photon ratio. The primordial abundance of 7 Li, as inferred from lithium 
abundances in Pop II stars and the analysis of stellar lithium depletion presented in this paper, is shown as 
a horizontal cross-hatched band. Vertical cross-hatched bands represent the allowed ranges of r\ based on 
the primordial abundance of 7 Li. 
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